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ABSTRACT 

We perform general relativistic ray-tracing calculations of the transfer of po- 
larized synchrotron radiation through the relativistic accretion flow in Sagittarius 
(Sgr) A*. Based on a two-temperature magneto-rotational-instability (MRI) in- 
duced accretion mode, the birefringence effects are treated self-consistently. By 
fitting the spectrum and polarization of Sgr A* from millimeter to near-infrared 
bands, we are able to not only constrain the basic parameters related to the MRI 
and the electron heating rate, but also limit the orientation of the accretion torus. 
These constraints lead to unique polarimetric images, which may be compared 
with future millimeter and sub-millimeter VLBI observations. In combination 
with general relativistic MHD simulations, the model has the potential to test 
the MRI with observations of Sgr A*. 

Subject headings: black hole physics — Galaxy: center — plasmas — polariza- 
tion — radiative transfer — sub-millimeter 



1. Introduction 

The compact radio source Sagittarius (Sgr) A* at the Galactic center is asso ciated with 



a ~ 4 x 1O 6 M supermassive black hole (jSchodel et al.ll2002l ; iGhez et al.ll2005l ) and is one 



of the best astrophysical sources for studying physics of black hole accretion. Since the 
discovery of active galactic nuclei and galactic X-ray binaries, it has been suggested that 
black hole accretion be responsible for the high luminosity and radiation efficiency of these 
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sources. To release the gravitational energy of materials falling into black holes, the angular 
momentum must be transported to large radii in the accretion processes. The nature of 
the viscosity resp onsible for the angular mom entum transport and energy dissipation is not 
well understood. IShakura k, Sunyaevl (119731 ) proposed that the viscous stress is induced by 
turbulence and is proportional to the pressure of the accretion flow. 

MHD simulations have made the magneto-rotational-instability (M RI) the primary 
mechanism for generation of turbulence and viscous stress in accretion flows (IBalbus fc Hawley 
199ll . Il998l ). Although there are indications that the viscous stress is proportional to the 
magnetic field energy density, the ratio of the magnetic field energy density to the gas pres- 



sure d epends on the init ial magnetic field configurations of these simulations (jPessah et al. 



2006 


)■ 


Kine et al. 


(2007) 



MRI may not be able to generate high enough viscosity to account for the observations. Nev- 
ertheless, the MRI is believed to play dominant roles in some black hole accretion systems. 
It has been studied extensively through MHD simulations and is awaiting observational test. 



Sgr A* plays a crucial role in this aspect. Since discovered in 1974 (IBalick fc Brown 



1974 ). the source has been obser ved extensively due to its association with the supermassive 
black hole iMelia fc Falckell200lh . Radio imaging reveals a co mpact radio source with its in- 



trinsi c size marginally uncovered in the millimeter (mm) band (IBower et al.ll2004l ; IShen et al. 



20051 ). The centimeter (cm) spectrum can be fitted with a power-law f unction, and there is 
clear evidence for a mm and sub-mm spectral bump (IFalcke et al.lll998l ). Radio (< 43 GHz) 
emission of Sgr A* has unusual polarization characteristics with no detection (< 0.4%) 
of linear polari zation (LP) and w eak circular polarization (CP) of up to ~ 1.0% in the 
quiescent- state ( IBower et al.l 120021). Weak LP (< 1%) wa s detected recent l y at 43 and 22 
GHz during radio outbursts ( Yusef-Zadeh et al. 2007al lbl). iMacquart et al. ( 2006 ) reported 
a mean LP of 2.1 ± 0.1% at 3.5 mm. In the sub-mm band , high LP of ~ 10% has been 
routinely measu red (Aitken et al. 2000 ; Marrone et al.|[2007h without any firm indication of 
detectable CP (IMarrone et al.ll2006l ). These observations suggest that Sgr A* is powered 
by accretion of the supermassive black hole with the mm and sub-mm emission originating 
within 10 R s , where R s ^ 1.2 x 10 12 [M/(4 x 10 6 M Q )] cm is the Schwarzschild radius for 
a non-spinning black hole w ith mass M. The high LP implies a very low accretion rate 
( Quataert fc Gruzinovl 2000 ). Based on the MRI, a small (< 10 R$) accretion torus can 



account for the spec tr um and po l arizat ion of the sub-mm emission (IMelia et al.l |2000| . 12001 



Bromley et al.l l200ll ). iLiu et aL ( 2007 ). hereafter L07, recently generalize the model to a 
two-temperature MRI driven accretion flow, which naturally explains the 90° flip between 
the electri c vector position angles (EVPA) of th e polarized sub-mm a nd near-infrared (NIR) 
emissions (jEckart et al.ll2006l ; iMeyer et al.ll2007l ; IMarrone et al.ll2007l ). Besides the basic pa- 
rameters related to the MRI, the model introduces only one more parameter to characterized 
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the not-well-understood electron heating processes. 

However, all models proposed so far have ignored the birefringence effects of synchrotron 
radiation propagating through the relativistic accretion flow in Sgr A* and have treated the 
general relativistic (GR) effects with simplified disk structures. The Faraday rotation of the 
birefringence effects can depolarize low frequency emission, and the GR light bending effect 
can complicate the transfer of polarized emission near the black hole. Based on the model 
proposed by L07, in this Letter we present a more complete and self-consistent treatment of 
radiation from the relativistic accretion flow. The ratio of the magnetic field to gas pressure, 
emission spectrum, and polarization are very sensitive to the onset frequency of LP from the 
cm to mm bands, electron heating rate, and disk inclination angle, respectively. A fitting to 
these observations leads to unique images of fractional LP and CP. The combined effects of 
GR and birefringence produce distinct LP and CP features between the sub-mm and NIR 
bands. These predictions can be tested with future polarimetric VLB1 observations. More 
importantly, in combination with MHD simulations, observations of Sgr A* can be used to 
test the MRL In §[2j we briefly discuss our treatment of the transfer of polarized synchrotron 
radiation near a Schwarzschild black hole. Our main results are presented in § [3j We draw 
conclusions and propose future work in § SJ 



2. Transfer of Synchrotron Radiation 



Following lLandi DeglTnnocenti fc Landoml (120041 ). hereafter LL04, we defined four vec- 
tors p, e, if, and p to treat the transfer of polarized emission in a magnetized plasma. Here, 
p represents three normalized Stokes parameters (Q, U, V)/I, and e,ff, and p correspond to 
three normalized emission coefficients (€q,€jj,€v) = (sq,£u,£v)/I, absorption coefficients 
(VQj Vu, Vv), and Faraday conversion and rotation (pq, pu, pv), respectively. For the to- 
tal specific intensity /, the corresponding emission and average absorption coefficients are 
indicated by ej = ejl and rj I , respectively. Then we have 

-(Vi + ff -p- e T )I, (1) 



ds 
dp 



—fj + (ff ■ p)p + p x p + e 
where s is the distance along a light ray. 



eip. 



(2) 



Synchrotron radiation is elliptically polarized with a dominant extraordinary emission 
component, whose electric field is perpendicular to both the magnetic field vector and the 



(Bekefi 


1966; 


Melrose 


1971) 



(4 + 4) 1/2 » ey 
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components as the axes, the U (or Q) components of e, ff, and p vanish, and pq (or pr/) 
and py correspond to the Farad ay conversion and rotation coefficients, respe ctively (Melrose 



1997; 



Kennett k Melrose 1998. hereafter M97;KM98). We use formulas of Melros 



e! Jl97lh 



to obtain e/ and e. Then rjj and ff can be derived from the Kirchhoff 's law. The polarization 
of the natural wave modes of a magnetized plasma determines p. Because there are no 
birefringence effects for the natural wave modes, p needs to be parallel to e (M97; KM98; 
LL04). For the thermal plasmas we are interested in here, M97 has given expressions for the 
Faraday rotation coefficients in cold and relativistic limits. We extrapolate these results py = 
(e 3 nB cos6 l J B/7r7 c mgC 2 z/ 2 )[7 c r 1 (l-ln7 c /(27 c )) + ln7 c /(27 c )], where e, m e , j c m e c 2 /k B , n, k B , c, 
B, 6b, and v are the electron charge, mass, temperature, density, Boltzmann constant, speed 
of light, magnetic field, angle between the magnetic field and wave propagation direction, 
and wave frequency, respectively. The Faraday conversion coefficient is then given by Pv{ e u J r 

4) 1/2 



ley. 



We simulate the pseudo-Newtonian disk structure in a Schwarzschild metric with all 
quantities treated as measured by a locally static observer. Because the emission frequencies 
in the mm to NIR bands are much higher than the cyclotron and plasma frequencies of the 
the rmal electrons in the ac c retion flow, the plasma effects on photon propagation discussed 
by iBroderick fc Blandfordl (120031 ) can be neglected and all photons travel along the null 
geodesies. To perform the abov e radiation transfer of polarized synchrotron emission with 
the ray-tracing code discussed in lHuang et al.l (120071 ). one only nee ds to take into accou nt the 
rotation of EVPA along a null geodesic due to the space curvature (IBromley et al.ll200ll ). The 
effects of magnetic field and G R have been taken care of in the radiation transfer equations 
(iBroderick fc Blandfordl I2004T ). 

MHD simulations of the MRI show that the magnetic field is dominated by its azimuthal 
component due to the shearing of the Keplerian accretion flow. Previous modeling of the 
sub-mm to NIR po larizations of Sgr A* assumed that the m agnetic field does not have other 
components, see in lMelia et al.l (120001 ) ; lBromley et al.l (l200ll );L07. Here we consider two more 
realistic configurations: one with the magnetic and velocity field parallel to each other in the 
upper half of the disk and anti-parallel in the lower half, and the other with the magnetic 
field direction just reversed. The averaged Stokes parameters of these two configurations 
are used to fit observations. Sgr A* is highly variable with variation timescales decreasing 
with obser ving frequency. In m m and sub-mm bands, both flux and LP vary on hourly 
timescales ( iMarrone et al.l 120061 ) . while the accretion timescale near the black hole is a few 
days. Shorter timescale variations are likely associated with turbulent fluctuations in the 
accretion flow, which may be addressed with MHD simulations. The results of our stationary 
accretion model should be compared with measurements averaged over days. 
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3. Spectral Fitting and Polarimetric Images 



The thick solid lines in Figfjja) and (b) represent the best fit to the spectrum and 
polarization data of Sgr A *. The data poin ts are summar ized in L07, except the two CP 
measurements at 112 GHz faower et alil2002h and 340 GHz iMarrone et al-lbood ). We fit the 
data aro und mm to sub-mm bu mp, including the EVPAs at 86 GHz averaged over five epochs 
in 2004 (IMacquart et al.l 120061 ). and at 2 30 GHz that rema ined unchanged during observa- 
tions from 2002 October to 2004 January (IBower et al.ll2005l ). In our fitti ng, neither the NIR 
measu rement ( Eckart et al.lboo"^ ) nor the detected LP at 22 and 43 GHz ( Yusef-Zadeh et al. 



2007cM is used simply because these are related to flaring activity. However, the EVPA at 
2.2/im remained unchanged from 2004 to 2006 and may be explai ned by a temporary disk 
with toroidal magnetic field (jEckart et al.ll2006l ; iMeyer et al.l 120071 ). The corresponding five 
model parameters are the ratio of viscous stress to magnetic field energy density % = 0.7, 
the ratio of magnetic field energy density to gas pressure (3 P = 0.4, the electron heating rate 
indicated by the dimensionless parameter C\ = 0.44, the accretion rate M = 0.95 x 10 _8 M Q 
yr _1 , and the inclination angle of the accretion flow i = 40°. The most prominent GR effect 
is the amplification of the emission area due to the light bending near the black hole (see 
Fig 12]). At higher frequencies the disk becomes optically thin and lower frequency emission 
is produced at relatively larger radii, so the flux density enhancement by the GR effects is 
significant near the spectral peak and vanishes in the low and high frequency limits as can 
be seen by comparing with the thin solid line of the pseudo-Newtonian spectrum. 

The GR effects also change the LP. The thin solid line in Figfjjb) shows the corre- 
sponding pseudo-Newtonian result. The LP is significantly suppressed (< 1%) by bire- 
fringence effects at low frequencies. Another radio component introduced by L07, which 
i s presumably associated with an outflo w, may fit the cm CP and flux density spectra 



outbursts ( Yusef-Zadeh et al 



2007a 



(IBeckert fc Falckeil2 002; Zh ao et al.ll2004f) and explain the polarization observed during radio 



b|). The dotted line shows the GR calculation without 
the birefringence effects, i.e. p = 0. We see that there is a weak level of LP (~ 2 — 3%) at low 
frequencies. The birefringence effects reduce the LP in mm and sub-mm bands by up to 50%. 
At very high frequencies, the birefringence effects are weak and the results approach that of 
the fiducial model. As shown by L07, the flux density is very sensitive to the electron heat- 
ing rate C\, which determines the electron temperature profile. However, the fitting to the 
emission spectrum doesn't give a unique set of model parameters. With the additional LP 
measurements, the best fit model parameters can be fixed. This is because the synchrotron 
emissivity is proportional to B 2 n while the Faraday rotation measure is proportional to Bn, 
for a B 2 n constrained by the emission spectrum, the Faraday rotation increases with the 
decrease of B, which is determined by (3 P . Therefore the onset frequency of LP increases 
with the decrease of /3 p . The two magnetic field configurations give almost identical LP but 
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different EVPA and CP. So we plot only one LP [dashed line in FigfjTb)] but two sets of 
EVPA and CP [dashed lines in FigfjTc) and (d)]. The average of Stokes parameters over 
the two configurations reduces the LP. One may use the difference between the average and 
individual configurations to estimate the variability of the corresponding quantities. 

The most interesting prediction of our model is oscillations of LP and CP with frequency 
~ 1 THz (see Fig{T](b) and (d)). This is due to the combined effects of GR and birefringence. 
Since p oc v~ 2 , the birefringence effects are negligible > 1 THz. The modulation of LP and 
CP by the birefringence is amplified by the GR effect ~ 1 THz, where J \p\ds is on the 
order of 2n. The LP oscillates with the amplitude and period decreasing due to rapid 
increase of the Faraday rotation angle as frequency decreases, in contrast to the pseudo- 
Newtonian and no-birefringence GR cases. This effect also gives rise to high CP of > 5% 
for each individual magnetic field configurations. Future polarization observations should 
be able to easily test these predictions. The model predicted intrinsic EVPA flips by ~ 90° 
three times from cm to NIR bands. The flip at ~ 50 GHz is related to the fact that the 
red-shifted side of the Keplerian flow becomes transparent and strongly polarized in the 
direction perpendicular to thedisk axis. At ~ 100 GHz, the near and far sides also become 
transparent and dominate the LP. The EVPA flips to be roughly aligned with the disk axis. 
At ~ 3 THz, the polarized emission is dominated by the blue-shifted side and the EVPA 
flips again. Without introducing an external Faraday rotation measure, the range limited 
by the two configurations fits the EVPA measurements well. The corresponding position 
angle of the disk axis is ~ 120° ± 30 °, which is different from that inferred from the NIR 



flare observations (IMeyer et al.l 120071 ) and the pseudo-Newtonian mo del (L07). The fitting 



to the EVPA can b e improved with an external rotation measure (IMacquart et al.l 12006 



Marrone et al.l 120071 ). which also brings the EVPA of the disk axis in agreement with the 
pseudo-Newtonian model. Similar to the pseudo-Newtonian case, the LP and EVPA flips 
are very sensitive to the inclination angle of the disk. In practice, cases with inclination 
angles ranging from 30° to 50° can fit the spectrum and LP within error bars. The best fit 
EVPA of 40° is 10° smaller than the pseudo-Newtonian result. Nearly edge-on (i > 60°) 
cases overestimate the LP in the sub-mm band and fail to produce a 90° EVPA flip from 
sub-mm to NIR bands because the near side of the torus dominates the emission. 

FigfS] shows the intrinsic polarimetric images of the fiducial model at 230 GHz with 
the disk axis in the vertical direction. Panels (a), (b), (c), and (d) are the images of the 
total specific intensity J, fractional LP (^/Q 2 + U 2 /I), fractional CP (V/I), and the residual 
(A0 res ) of the Faraday rotation angle J (pq + pfjY^ds divided by 2tt, respectively. Results 
shown in Fig{T] are obtained by integrating the Stokes parameters over the disk to compare 
with spatially unresolved observations. Doppler boosting enhances / on the left hand side 
of Figfj^a). The white dashes in Figj2](b) represent the intrinsic EVPA of each ray. In 
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Figj2^c) positive (red) and negative (black) values of CP correspond to right and left hand 
polarization, respectivehJE The strips in Figf5](b) and (c) are clearly correlated with the 
Faraday conversion shown in Figfj^d). In our case, the Faraday conversion coefficient is 
always much greater than the Faraday rotation coefficient, and its integrations along the 
rays are much larger than 2tt at 230 GHz, the Stokes V therefore changes its sign over 
one period of EVPA rotation (KM98; LL04). The strips will be less distinct in a realistic 
turbulent flow. The observed images will alway s be smoothed by finite spatial resolutions and 
interstellar scattering effects (IShen et al.ll2005l ). However, rough contours may be resolved by 
future polarization-sensitive VLBI observations since the scattering effect is less important 
in sub-mm band. The black hole shadow may be detected in the / image, as evidence for the 
GR theory. We also expect high levels of LP with the EVPA mostly in the radial direction, 
as evidence for toroidal component dominated magnetic fields. 



Discussion and Conclusion 



Treme ndous efforts have been put to image shadows caused by the strong gravity of 
black holes (jKrichbaum et al.ll2006l ; IShen et al.ll2005l ; iBower et al.ll2004l ) . Sgr A* has been the 
best candidate for such studies. The polarization measurements of Sgr A* have put strong 
constraints on the physical processes near the black hole. Based on the two-temperature 
MRI driven Keplerian accretion model, which has only five parameters, we calculate the 
synchrotron emission from the accretion flow with both the GR and birefringence effects 
taken into account self-consistently. The model can reproduce the spectrum and polarization 
of Sgr A* in mm and shorter wavelength bands and further predicts distinct LP and CP 
at ~ 1 THz. Polarimetric images of the relativistic accretion flow of the best fit model are 
made, which can be verified with future polarization-sensitive VLBI observations. 

The MRI has become the leading mechanism for the generation of viscous stress required 
to drive the accretion flow of plasmas into black holes. The theory has been well established 
through both theoretical analyses and numerical investigations. The model presented here 
makes it possible for the first ti me to put the theory into observational tests. If the scaling 
laws for the MRI discovered by iPessah et al.l (120061 ) are applicable to the accretion flow in 
Sgr A*, the fiducial model parameter (3 P = OA requires the Alfven velocity induced by the 
large scale vertical magnetic field be ~ 20 times smaller than the sound velocity. With 



6 The convention used for CP (V/I) is consistent with the IEEE definition which has been adopted by 
the IAU and used by many observations, i.e., right hand polarization corresponds to the electric vector of 
an approaching wave rotating counter-clockwise. 
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this lar ge scale magnetic field configuration, one can carry out GR MHD simulations for 
Sgr A* (INobel et al.ll2006l ). The outflows from these simulations should not only reproduce 
low frequency radio flux, CP, and source size measurements, but also induce a Faraday 
rotation measure external to the disk studied here, which can improve the fit to the observed 
EVPAs in mm and sub-mm bands. NIR and X-ray flares and the high variability of this 
source, whic h cannot be studied with the current model, may also be addressed with such 



simul ations (IBaganoff et al.l l200ll ; iGillessen et al.l 120061 ; iHornstein et al.l 120071 ; iMeyer et al. 
2007h . 
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Fig. 1. — Plots of the spectrum, LP, EV PA and CP of Sg r A* in panels (a),(b),(c), and (d), 
respectively. All the data points are from lLiu et al.l (120071 ). except t he two CP measurem ents 
at 112 and 340 GHz (open triangles in panel (d)), which come from lBower et al.l (120021 ) and 



Marrone et al.l (120061 ). respectively. Although the NIR data are for flares, the optically thin 



emission is always dominated by emission from the blue-shifted side and the LP therefore 
can be compared with the model. The thick solid lines in panels (a) and (b) correspond to 
spectrum and LP of the fiducial model and the thin solid lines represent the corresponding 
pseudo-Newtonian results. The GR effects increase the flux density by a factor of ~ 2 near 
the spectral peak. The dashed lines in panels (b),(c), and (d) show LP, EVPA, and CP of 
two magnetic field configurations, respectively. For LP, the two configurations give almost 
identical results so that only one is drawn. The dotted line in panel (b) shows the result 
with no birefringence effect, i.e. p = . We predict strong CP variations near 1 THz and 
three 90° flips in the EVPA of the LP (see text for details). 
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Fig. 2. — Polarimetric images of the accretion flow in Sgr A* at 230 GHz. The length unit is 
Rs and the disk axis is in the vertical direction. Pane l (a) shows the total i ntensity map. The 
left h and side is brightened by the Doppler boosting (IBromley et al.ll200ll ; iBroderick fc Loeb 
20061 ). Panel (b) gives the map of the fractional LP with the EVPA and magnitude indicated 
by the white dashes and the color wedge, respectively. Panel (c) shows the fractional CP 
with the positive sign for right hand polarization (IEEE convention, see text for details). 
Panel (d) gives the residual of the integration of the Faraday conversion coefficient along the 
rays divided by 2tc, with which the strips in the LP and CP images are clearly correlated. 



